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We show that there are regions of parameter space in muhi-scalar doublet models 
where, in the first few hundred inverse femtobarns of data, the new charged and 
neutral scalars are not directly observable at the LHC and yet the Higgs decay 
rate to 66 is changed significantly from its standard model value. For a light Higgs 
with a mass less than 140 GeV, this can cause a large change in the number of 
two photon and r+r^ Higgs decay events expected at the LHC compared to the 
minimal standard model. In the models we consider, the principle of minimal fiavor 
violation is used to suppress fiavor changing neutral currents. This paper emphasizes 
the importance of measuring the properties of the Higgs boson at the LHC; for a 
range of parameters the model considered has new physics at the TeV scale that is 
invisible, in the first few hundred inverse femtobarns of integrated luminosity at the 
LHC, except indirectly through the measurement of Higgs boson properties. 
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I. INTRODUCTION 



Experiments at the LHC will directly probe physics at the weak scale. Most 
physicists believe that there is new physics at this energy scale, beyond what is in 
the minimal standard model (SM). This belief is motivated to a large extent by the 
hierarchy puzzle and by the fact that the scalar sector of the standard model has yet 
to be directly probed by experiment. The single doublet in the SM is the simplest 
example of a scalar sector but many extensions have been studied. Amongst the most 
widely considered are two Higgs doublet models ^. A problem that immediately arises 
in such models is the possibility of flavor changing neutral current (FCNC) effects 
that are unacceptably large. In particular, when the SM fermion fields couple to 
both doublets, and the couplings are arbitrary, FCNC effects are possible at tree 
level. Glashow and Weinberg gave a simple prescription for how to avoid such effects 
through imposing a discrete symmetry [2!|. One can also suppress FCNC effects by 



adopting an ansatz suppressing the coupling of the new doublet 
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In this 
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paper, we use the principle of minimal flavor violation (MFV) 
which causes tree level FCNC to vanish (or at least be suppressed by small mixing 
angles) in multi doublet models in a natural way. 

We are interested in the possible effects of new physics on the properties of the 
Higgs boson and we characterize the impact of new physics on the Higgs through 



an operator analysis 



is Q, 
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IgI . 17, Q, Q • We assume that the new physics 
mass scale M is much larger than the Higgs boson mass and add higher dimen- 
sion operators that are invariant under the S'f/(3) x SU{2) x U{1) symmetry of the 
standard model. Since the operators give small corrections to the SM, one does not 
expect them to influence standard model processes that are unsuppressed. For exam- 
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See m for a review. 
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pie, the Higgs coupling to two W bosons is an unsuppressed tree level coupling and 
new physics contributions to it should be negligible. However the dominant Higgs 
production mechanism through gluon fusion, g g ^ h occurs at leading order in per- 
turbation theory through a top quark loop. Hence, in the SM it is suppressed and 
new physics can easily compete with the standard model contribution 13|. Similar 
remarks hold for the /i — 77 decay amplitude. Some of the tree level couplings of 
the standard model Higgs are also very small. For example, the Higgs to r Yukawa 
coupling is of order^ mT-/v ~ 0.75 x 10^^ and the Higgs to b-quark Yukawa coupling is 
of order nib/v ~ 2 x 10^^. New physics characterized by higher dimension operators 
can also compete with the standard model in the h rf \v\ and h ^ bb decay 
amplitudes. 

A light Higgs with a mass less than 140 GeV is likely to be detected first through 
its decay to two photons despite the fact that the branching ratio for this process 
is quite small, i.e. of order 10~^. Early detection through its decay to r+r^, which 
has a branching ratio around 10^^, may also be possible. The dominant decay mode 
is to 6 6 pairs. However this decay mode is much harder to observe because of large 
theoretical uncertainties^ on the cross section of the irreducible SM background tibb. 
An integrated luminosity of order 500 fb~^ may be required to observe the standard 
model Higgs in the b b channel. If the rate for h ^ bb decay is changed by a factor / 
from its standard model value (but other properties of the Higgs are left unaltered) 
then the number of /i — 77 or /i — r^r^ decay events observed at the LHC is 



^ Here v ~ 250 GeV is the vacuum expectation value that spontaneously breaks the weak gauge 
group down to the electromagnetic gauge group 



20| for a recent review on production and detection of the Higgs and 2l| for a recent study 
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on the tibb SM background. 
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changed by a factor ^ where 

£ = = (I) 

l + (/-l)Br(/i^66)sM' 

and Br(/i &&)sm is the SM branching fraction of h ^ bb. Even though h bb 
decay will not be directly observable until there are many years of LHC data, it's 
rate is of crucial importance for all measurable properties of the low mass Higgs. 

In this paper, we consider multi-doublet scalar models.^ For simplicity we restrict 
our attention to two scalar doublets H and 5*. Here H denotes the usual SM Higgs 
doublet with a mass less than 140 GeV. 5* is a new scalar doublet with a mass 
M ~ 1 TeV and has the same quantum numbers as the SM Higgs doublet H. We 
demonstrate that there are regions in parameter space for which the new doublet 5* 
will be invisible at the LHC, at least in the first few hundred inverse femtobarns of 
data. The largest observable effect of this new doublet S will be order one shifts in 
the h ^ bb rate which in turn will affect the branching ratios for all light Higgs decay 
channels. It is straightforward to generalize our results to a scenario with more than 
two scalar doublets. 



* New physics in the form of a second scalar doublet with an unbroken Z2 symmetry can also 
provide a component of dark matter and the effect of such a second doublet on the SM Higgs 
was recently examined in 22 1 . 
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II. THE TWO DOUBLET MODEL 

The scalar potential for the model we consider is given by ^ 

V{H, S) = ^ (^H^H - y) ' + M'S^S + ^ (S^ sf - [g, [S^H) [H^H) + h.c] 

+02 [S^S) [WH) + [g'^ [S^H) (S^H) + h.c] + g'^ (S^H) (H^S) 

+ [g,{S^S){S^H) +h.c.]. (2) 

We assume the phase of S is adjusted so that gi is real. S appears with a positive 
mass term and acquires a vacuum expectation value only through it's coupling to H, 
which undergoes the usual electroweak symmetry breaking. Since M is much greater 
than the weak scale v, the neutral component gets a vacuum expectation value 
that is much smaller than v 

(50) ^ -^1^ « V. (3) 

In addition to the SM Yukawa couplings of the doublet if, the doublet S has the 
following Yukawa couplings to the quarks, 

A£y = -riDdndDS'^QL - VuUnguSeQL + h.c. (4) 

We make use of the principal of minimal flavor violation. This results in small loop 
level FCNC through the appearance of Yukawa coupling matrices gr) and gu in the 
loops. We assume that possible multiple insertions of the Yukawa matrices in Eq. (jl]) 
are suppressed. In this model, the physical quark masses are the result of the sum of 
contributions from the coupling of quarks to H and S. Thus in the mass eigenstate 
basis, the Yukawa matrices gu,D do not satisfy the usual relation to physical quark 



^ We thank Lisa Randall for pointing out the typo in the sign of gi in Eqn.(2) in the previous 
version [23 1. 
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masses. In the SM, g\j j^, = \/2mi/v. In the down quark sector, in the mass eigenstate 
basis, the couplings of the heavy scalar doublet are 



^Cy = -7]DV2dR — V^ULS' - 7]DV2dR — dLS'^ + h.c. (5) 
V V 



Here V is the CKM matrix, go = V^m-d/v and nid is related to the physical down 
quark mass by 

irid = . 6) 

We assume that the constant rju in Eq. (jlj) is very small so that the S coupling 
to the up-type quarks can be neglected. When rju <ti I the production of S via it's 
coupling to the top quarks is suppressed. 

On the other hand, we take rio to be large, r^/j ~ 10 and for simplicity we choose 
it to be real. Since the down type quark Yukawa couplings in the matrix are very 
small, the effective coupling rjogD is still perturbative. The choice of rjD ^ 1 makes 
the coupling of b quarks stronger to S compared to H resulting in a large shift in 
h ^ bb once S is integrated out. Thus, with rju ^ I, rju 1, and M t> an almost 
invisible S can be produced at the LHC and will leave it's footprint through large 
shifts in the h —>■ bb rate. 



III. FCNC CONSTRAINTS 

Even though we imposed MFV eliminating the possibility of tree level FCNC, in 
our two doublet model, there are at one loop corrections to standard model FCNC 
processes. In particular, we want to check that the choice of Irjol ~ 10 ^ 1 is 
consistent with constraints from FCNC. Consider the weak radiative b decay b ^ s'y 



with two doublets 24] • Calculating the Feynman diagrams in Fig.([T]) we find that 






FIG. 1: The one loop contribution to 6 ^ 57 due to the doublet S. 



charged S exchange induces at one loop the effective Hamiltonian, 



(7) 



where e < is the electron charge. 

For inclusive decay, this does not interfere with the standard model contribution 

I; 

from O7 [25] since the strange quarks have opposite chirality. Hence, we find 



F(6 ^ 57) 



1 + 



r/l) Tils rrib mf 



(8) 



Taking r]D = 10, c/i = 1, M = 1 TeV and \C7{mb)\ ~ 0.3 in Eq. §^ gives AF(6 ^ 
S7)/F(6 s7)sM ~ 10^^, which is much too small to be observed^. For exclusive 
decays there can be interference between the standard model contribution and the 



contribution of the effective Hamiltonian in Eg . (jT l 
cloud our ability to constrain the new physics 



but there hadronic uncertainties 



33| 



There are contributions from the interactions of the new doublet S to the Wil- 
son coefficient C7 that are proportional to rjoVu and are not suppressed by rhs/mb- 



26 



27| for the latest calculation of BR(i3 Xg 7) at NNLO in QCD and its comparison to 



See 



the results from Babar 



28l |29| and Belle 30| as averaged by the Heavy Flavor Averaging Group 



3l|. The remaining uncertainties in theory and experiment preclude a exclusion of our model 



based on 6 — *■ 57 constraints for the parameter space of interest. 
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However, because we have focused on the region of parameter space where rju is very 
small these have been omitted. Our assumption of very small rju greatly diminishes 
the constraint that b ^ s'~f places on the model. 



IV. EFFECTS ON LIGHT HIGGS DECAYS 



So far we have described the general features of the two doublet model we are 
considering and demonstrated the compatibility of a large |?7£)| ~ 10 with 6—^57 
constraints. Next we investigate how this simple extension of the SM affects light 
Higgs decay to quarks by integrating the heavy doublet S out of the theory to induce 
the effective operator 



= -VD{H^H)j^dRgDH^QL + h-c 



Including the effects of this operator and using Eqs. 
h ^ bb rate is modified relative to the SM as 



and 



r{h bb) 
r{h^ bb) 



SM 



1 + 3v^gir]D/2M^ 



(9) 

we find that the 
(10) 



1 + v^gir]n/2M^ 

Note we have included terms suppressed by more powers of f ^/M^ in Eq. (fTOj) which 
are accompanied by the large factor r/^. However, we can still consistently ignore 
the effects of dimension 8 operators contributing to h ^ bb since their contributions 
start at order ri£, v'^/M^ <^ 1. 

At M = 1 TeV, the parameter choices of gi = 0.5, rjo = 10 or gi = l,riD = 5 give 
the rate for h ^ bb which is 1.6 times it's SM value. An even more dramatic effect is 
seen for the parameter choices of gi = —2, rjD = 20 and gi = 1, r//) = —10 which give 
a rate that is 121 and 0.008 times the SM value respectively. Thus, the presence of 
an additional TeV scale scalar doublet with a coupling to b quarks about ten times 
the SM value can have dramatic changes in the decay width and branching fractions 
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of a light Higgs. For example, with gi = — 2,?7£) = 20 the branching ratio for the 
experimentally promising modes of — > 77 and h t^t^ will be down by a factor 
of ^ ~ 1/80. Such a scenario would make detection of the light Higgs very difficult. 
On the other hand for gi = l,rir) = —10 in which case the h ^ bb rate is 0.008 
it's SM value. In this case, the branching fractions for /i — 77 and h t^t^ will 
increase by the factor ^ ~ 3 for mh = 120 GeV. This would apply for Higgs searches 



at the Tevatron as well 3J]. 

One can also generahze the above analysis to the lepton sector and induce a 
corresponding effective operator 

£efr = -m{H^H)j^eR~g,H^LL + h.c, (11) 

which contributes to the decay of the Higgs to charged leptons h — > By 
choosing a large value for rji one can similarly induce order one shifts in the decay 
rate to h ^ £^£^ . Such order one shifts can be seen at the LHC in the experimentally 
promising channel of h t^t^ . The effect of the operator in Eg. (11 II) was recently 
studied in 17| where naturalness criteria were used to constrain the size of the Wilson 
coefficient. It was shown that order one shifts are indeed possible and compatible 
with experimental constraints. The branching ratio for h t^t^ can be influenced 
both by the effect of the operator in Eq. ( ITTl) on the rate for h — > r^r^ and by the 
effect of the operator in Eq. IQ on the total Higgs width. 

Order one shifts in the rate for h t'^t' can also affect the total width of the 
light Higgs since it's branching ratio is not negligible. For example, if rrih = 120 
GeV then the branching ratio in the SM for h t^t^ is about 7%. As order one 
corrections are possible to partial decay widths for both h bb and h t^t^ , the 



New physics in the form of a massive fourth generation neutrino or additional scalar singlets 
36| have also been shown to effect the possibility of the detection of the Higgs at LHC. 
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relative impact of the two decays on the total width can be changed dramatically. 
This can make the total decay width even more sensitive to the effects of S. However, 
for the sake of simplicity we will assume in this paper that rji is small so that effects 
on the width of the Higgs from the coupling of S to leptons are negligible. 

The number of observed h ^ '-fj events can also be effected by higher dimension 
operators that induce a direct coupling between h and 77 and between h and gg as 



discussed in 



13|. The latter effects the Higgs production rate by gluon fusion. New 



physics effects of this form are distinguishable from a change in the total width as 
the new physics effects on the total width will cancel in the ratios of the number of 
expected events for different Higgs production mechanisms and decay channels. 



V. PRODUCTION AND DECAY OF THE NEW SCALAR DOUBLET 

We now study the production and decay of the new scalar doublet S and discuss 
the possibility for it's observation at the LHC. The doublet contains new neutral 
and charged scalars with masses approximately equal to M. If M ~ 1 TeV the LHC 
has enough energy to produce these states. However, we will show that for a range 
of parameters their production rates are quite small and that the dominant decay 
channels have poor experimental signatures making them invisible at the LHC, at 
least for the first few hundred inverse femtobarns of data. 

The production of the charged S"^ is suppressed compared with the neutrals and 
the pseudoscalar Sj does not have a significant branching ratio to the most promising 
detection channels WW and ZZ. Hence we present in detail a discussion of the 
neutral scalar^. Expanding H and S about their vacuum expectation values we 



Similar conclusions hold for the charged scalar and neutral pseudoscalar. 
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write 

qO I -c-O 

- (H^) = ^, S'- {S') = (12) 

The fields and mix and the resulting mass eigenstate fields h and Sr are 
approximately given by^ 

hc^h' + Sn^S'^- ^^h'. (13) 

The production rate for Sr is very small. The dominant production mode is 
through bb ^ Sr. In this process the initial b and b each come mostly from coUinear 
gluon splitting and the remaining spectator b quarks have very low transverse mo- 
mentum to be observed in the final state. The large logarithms associated with 



collinear gluon splitting into light c 
by one or two orders of magnitude 



uark pairs leads to an enhancement of the rate 



37 
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391 ] over g g ^ bb Sr^^ where the scalar 



is radiated off one of the final state 5-quarks. This is because the final state 6-quark 
which radiates the scalar is far offshell before emission and thus the rate does not 
receive the enhancement of large logarithms associated with collinear gluon splitting. 
The cross-section for 6 6 ^ Sr at leading log takes the form 

a{bb^SR)^^c^'^(%) f ^b{x,^,)^^,^,), (14) 
3 s V ^ / JAiys X XS 

where 6(x, /i) and fe(x, /i) are the b quark and antiquark parton distribution functions 
respectively and s is the center of mass energy squared. The large logs from collinear 
gluon splitting are summed into the parton distribution functions by choosing ~ M. 
As seen from Eq. (1141) , the bb Sr cross section receives an additional enhancement 
by a factor of r/l, ~ 100 compared to the production of a SM Higgs with the same 
mass. This production cross section as a function of the mass M is shown in Fig. ([2]) 



^ We assume that the parameters in the scalar potential are real so there is rin t ,-''? juixing. 
^° This result is based on the NLO QCD calculation for tt + h production in [4d[4lU4a. |43|. 
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FIG. 2: The production cross section of 66 ^ Sr for the parameter choices of 7]£) = 10 
and gi = 0.5. The production cross section (for the same parameter choices) for gg ^ Spi 
is approximately 0.2 fb for M ~ 800 GeV and falls quickly with increasing M and thus 
is not shown. The curve was generated using CTEQ5 parton distribution functions [44]. 
The 6 and t masses were evaluated at /i = 1 TeV using leading log running. We chose 
^QCD = 0.1 GeV and used the initial values m.}, = 4.26 GeV determined from converting 
the result of the IS* fit to the b quark mass [45] and the value nit = 170 GeV from the 



as the solid black curve. This curve was generated for the choice of riD = 10 and 
Qi = 0.5. We see that at M = 1 TeV the bb Sr cross section is about 10 fb. Thus, 
for 100 fb~^ of data one can expect the production of about 1000 neutral scalars Sr 
from b b fusion. Note that this dominant production mechanism doesn't exist for the 
heavy charged scalars. 

The next largest production mode of Sr is through gluon fusion and is given by 



PDG [46i] 
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a direct modification of tlie SM cross section 

M2 



47| 



(^{99 'S'ij)LHC 



a2 



6Ans v"^ 



rrib 
Vd — I 
nib V 



Ami 



2M2 V4m2 



M2 



F[/i,M,s], (15) 



wfiere we liave used tlie functions 

F[/i,M,s] = 



—9{^,N9\ ,/^), 



(16) 



and \{y) for y > 1 which is given by 48j 
1 



\_ y 

2y^ 2y2 



iTT log [y/y + y^y~l) - log^ ( y/y + y^y - l) + 



21 



(17) 



Here g{x, /i) denotes the gluon parton distribution function. As seen in Eq. (fT5|) this 
production channel receives significant contributions from bottom and top loops. 
The bottom loop has a significant contribution due to the direct coupling of Sr 
which involves 77^) ~ 10 ^ 1. Although the direct coupling of S*^ to the top quark 
is negligible for r/^/ -C 1, the top loop still gives a significant contribution due to the 
mixing of with the Higgs With the same parameters used as in Fig. ([2]) at 
M = 1 TeV one can expect the production via gluon fusion of only about 8 neutral 
scalar s Sr for 100 fb~^ of data. 

Other production mechanisms are similarly small. For example Higgs production 
(via vector-boson fusion) in association with massless jets, qq qqh, where q = 
{u, d, s}, is dominated by the Higgs being radiated off a virtual W or Z boson. So 



a{qq^ qqSn) 



V2M2 



< 1. 



(18) 



a{qq -^qqh)sM 

The pattern of possible decays of the new scalars in the model depends on the 
mass splittings between the various states. To simplify our discussion of the spectrum 
we neglect the S*^ — mixing and assume that the coupling constants in the scalar 
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potential are real. Then there is no — Sj^ mixing and the mass spectrum is, 



^s± - M'^ + g2 v"^ 
ml 

4? 



m%^M^ + {g2-g'2 + g'il2)v\ (19) 



We focus on the region of parameter space where the lightest scalar is 5*^. For its 
decays it is important to include the effects of 5*^ — mixing. The most important 
decay modes of Sr have the partial rates 

-(^-"-)-iii(^)^ (-) 

2 2 

TiSR-^W+W~) (22) 
^ ' 167rM' ^ ' 

TiSn^ZZ) ~ ^i^, (23) 

r(S«-/.w.)-^. (26) 

A plot of the branching fractions of Sr as a function of the mass M is shown in 
Fig. ([3]) for (7i = 0.5, M = 1 TeV and rjj^ = 10. The dominant decay channels 
are Sr hh and Sr bb. The Sr — > bb channel is known to have a large SM 
background. As we will discuss later, even the S hh channel can be difficult to 
observe. The final states where the Sr decays to gauge bosons and at least one 
of the gauge bosons decays to electrons and/or muons have a cleaner experimental 
signature. Note that for the parameters used in Fig. ([3]) the total width of an Sr 
scalar of mass 1 TeV is only 3 GeV. For comparison, note that the width of a 
standard model Higgs with a mass of 1 TeV is about 700 GeV 



491]. This is because 
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Branching Fractions for Sr Decays 
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FIG. 3: Branching fractions for Sr decays as a function of the it's mass M. The solid black 
curve denotes the branching ratio for 5/j bb, the gray very-fong-dashed curve denotes 
Sr hh, the red short-dashed curve is for Sr ^ tt, the blue medium-dashed curve is for 
Sr. W^W~ , and the green long-dashed curve is for Sr Z^Z^. We have not shown 
the curve for Sr hhh in order to avoid too much clutter. These curves were generated 
with the parameter choices of r]D = 10 and gi = 0.5. 

of the small vacuum expectation value of the heavy doublet which suppresses the 
coupling of Sr to two gauge bosons. 

In Table. [T]we show the number of expected events for 100 fb~^ of data at the LHC 
when Sr decays to gauge bosons and M = 1 TeV. For example with gi = 0.5, M = 
1 TeV and 7]d = 10 we find that a{pp SRX)Br{SR W+W')Br{W+W- 
i^ui^ui) ~ 0.05 fb, where we have summed over I = e,n leading to about five 
events with 100 fb~^ of data. For these parameters, the heavy scalar Sr will not be 
detected at the LHC in the first few hundred femtobarns of integrated luminosity. 
In fact, within much of the region of parameter space where the coupling of the new 
S-doublet to charge 2/3-quarks is suppressed (i.e., rju very small) the heavy scalar 
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Decay Channel 


gi = 0.5, riD = 10 


gi = l,r]D = 5 


5i = -2,r/fl = 20 


Sr Z^Z^ 


0.23 


0.10 


8.0 


Sr Z^Z^ ^ i+£-uu 


1.4 


0.58 


47 


Sr w+w- e+ui-u 


4.6 


2.0 


160 


Sr^W+W- ^ {i+i^jj,e~iyjj) 


3.0 X 10 


12 


1.0 X 10^ 



TABLE I: Expected number of events for 100 fb^^ of data at the LHC in the experimentally 
favored decay modes of Sr for different choices of the parameters gi,r]D. The production 
cross section is the sum of the a{g g 5'/j)lhc and a{bb Sr)i,uc cross sections. We use 
£ to denote either an electron or muon (i.e, we have summed over I = e, ^) and j denotes 
a single jet. We have chosen the mass of Sr at M = 1 TeV and a center of mass energy of 
14 TeV. After realistic selection cuts the final number of accepted events will be lower. 



Decay Channel 


91 = 0.5, r]D = 10 


51 = l,??!? = 5 


gi = -2,r/B = 20 


Sr —> hh ^ 6677 


1.1 


0.46 


0.82 


Sr ^ hh ^ bbr^T^ 


34 


14 


26 



TABLE IL Expected number of events for 100 fb^ of data at the LHC in the experimen- 
tally favored channels when Sr first decays to a pair of light Higgses 5^. We have chosen 
the mass of Sr at M = 1 TeV and the Higgs mass at rrih = 120 GeV. 



degrees of freedom associated with the doublet S are difficult to detect at the LHC 
as shown in in the first two columns of Table. [B However, as seen in the third column 
of Table. [H even with rju small, there are regions of parameter space that are more 
promising for detection at LHC. For example, for gi = —2, M = 1 TeV and rjo = 20 
we find that a{pp SRX)Br{SR W+W-)Br{W+W- £+vl-Vi) ~ 1.6 fb 
and detection of the new heavy scalar Sr at the LHC with a few hundred inverse 
femtobarns of integrated luminosity is more likely. 
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In Table. [TTl "we show the number of expected events when Sr decays to a pair 
of hght Higgses and one of them decays to the experimentally favored 77 or r+r~ 
channels. As seen in the table, for the parameters chosen detection is unlikely. The 
number of events in the last column are suppressed because for these parameters the 
Higgs decay rate to bb is enhanced by a factor of about 100 which reduces the Higgs 
branching ratio to 77 and r+r^ by a similar factor. 

VI. CONCLUDING REMARKS 

We have demonstrated that for regions of parameter space in multi doublet models 
the states of the new doublets are impossible to directly detect at LHC, using the 
first few hundred inverse femtobarns of data, and yet the effect of the new doublet 
on the total width of the light Higgs is very significant. In the simple two doublet 
model we considered in detail, the promising 7 7 and h t^t^ signals at the 
LHC for detecting a light Higgs could be significantly enhanced or suppressed. This 
demonstration emphasizes the importance of determining the properties of the Higgs 
boson in the presence of new physics that is difficult to directly detect at LHC. 

We thank Marat Gataullin for many helpful comments. This work was supported 
in part by the DOE grant number DE-FG03-92ER40701 and DE-FG03-97ER40546. 
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